Does solar structure vary with solar magnetic activity? 
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ABSTRACT 

We present evidence that solar structure changes with changes in solar ac- 
tivity. We find that the adiabatic index, Fi, changes near the second helium 
ionization, i.e., at a depth of about O.QSR©. We believe that this change is a 
result of the change in the effective equation of state caused by magnetic fields. 
Inversions should be able to detect the changes in Fi if mode sets with reliable 
and precise high-degree modes are available. 

Subject headings: Sun: interior — Sun: activity — Sun: oscillations — Sun: 
helioseismology 

1. Introduction 

It is well known that frequencies of solar oscillations vary with time (Elsworth et al. 1990; 
Libbrecht & Woodard 1990; Howe, Komm & Hill 1999, etc.). These variations are also 
known to be correlated with solar activity. However, these frequency changes appear to 
result from variations in solar structure close to the surface, and no evidence has been 
found of structural changes deep in the solar interior (Basu 2002; Eff-Darwich et al. 2002; 
Vorontsov 2001; Monteiro et al. 2000, etc.). This is in contrast to solar dynamics, which 
shows changes correlated with solar activity (Basu & Antia 2003; Vorontsov et al. 2002; 
Howe et al. 2000; and references therein). Basu (2002) showed that the currently available 
mode sets do not allow us to make any statement about changes in solar structure above 
about O.95R0 through inversions; the lack of high degree modes makes the resolution of the 
inversions in these layers poor and errors in the inversion results large. 

There have been recent studies of solar active regions using local helioseismology tech- 
niques that suggest that active regions have lower sound speed compared with quiet regions 
just below the surface (Basu, et al. 2004; Kosovichev et al. 2000, 2001). It has also been 
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shown by Basu et al. (2004) that the adiabatic index, Fi of active regions is considerably dif- 
ferent from that of quiet regions and the amount of change increases with increasing strength 
of the active region. Figure 1 shows the adiabatic index for two of the active regions studied 
by Basu et al. (2004), and one can see that the magnitude of the depression in Fi at the 
second helium ionization zone decreases with increasing strength of the active region. 

This leads us to question whether similar changes occur globally in the Sun as solar 
activity increases — of course one would expect the changes to be much smaller than those 
seen in active regions because global activity levels are much smaller than in active regions. 
Although the He 11 ionization zone is too shallow for studying temporal changes directly by 
inverting the currently available data sets, there are indirect ways by which we can study 
this region. 

Any spherically symmetric localized, sharp feature or discontinuity in the Sun's internal 
structure leaves a definite signature on the solar p-mode frequencies. Gough (1990) showed 
that abrupt changes of this type contribute a characteristic oscillatory component to the 
frequencies ly^/ of those modes which penetrate below the localized perturbation. The am- 
plitude of the oscillations increases with increasing "severity" of the discontinuity, and the 
wavelength of the oscillation is essentially the acoustic depth of the sharp-feature. Solar 
modes encounter two such features, the base of the convection zone (henceforth CZ) and 
the He 11 ionization zone. The transition of the temperature gradient from the adiabatic 
to radiative values at the CZ base gives rise to the oscillatory signal in frequencies of all 
modes which penetrate below the CZ base; the depression in the adiabatic index Fi in the 
He 11 ionization zone causes the second signal. The two oscillatory signals have very different 
wavelengths and hence can be decoupled. Not all modes see both features, only low degree 
modes {£ < 25) see the CZ base, higher degree modes only see the ionization zone, and the 
very high degree modes see neither. This signal has been used previously to study the CZ 
base (Monteiro, Christensen-Dalsgaard & Thompson 1994; Basu, Antia & Narasimha 1994), 
and also study changes in that region (Monteiro et al. 2000). 

In this work we use the oscillatory signal from the depression in Fi to study whether or 
not there are changes in Fi in that region that are correlated with solar activity. 

2. Technique 

The amplitude of the signals from the CZ base as well as from the He II ionization zone 
are small, hence we first amplify them by taking the fourth differences of the frequencies: 
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Taking the fourth differences amphfies the signals enough to be able to isolate them, keeping 
the errors at a manageable level (see Basu, Antia & Narasimha 1994 for a discussion). 

The fourth differences can be fitted to the functional form given by Basu (1997): 

S'^u = [ai + a2iy + (as + a^L)/^'^] + 

f as + ^ + ^ + ^) sm{2iy^Ta + ^a) + 

\ '^m m m / 

(&i + ^ + ^ + ^) sm(2u^n + ^b), (2) 

\ / 

where 

L = i{i + 1), z/^ = z/ - 7aL/2rat/, and, = - jbL / 2TbU . (3) 

The coefficients 01-04 define an overall smooth term, 05-O8, Ta and ipa define the os- 
cillatory contribution due to the He II ionization zone and the remaining terms define the 
oscillatory contribution from the CZ base. The frequencies of the two components, and 
are approximately the acoustic depths of the He II ionization zone and the CZ base respec- 
tively, but they also include a contribution from the frequency dependent part of the phases 
ijja and i/jf, which is not taken into account explicitly. The coefficients in this expression are 
determined by a least squares fit to the fourth differences. 

We use oscillation frequencies determined by the Michelson Doppler Imager (MDI) on 
board the SOHO spacecraft, as well as data from the Global Oscillations Network Group 
(GONG). We have used 38 data sets from MDI (Schou 1999), each covering a period of 72 
days starting from 1996 May 1 and ending on 2004 March 19. We use 28 data sets from 
GONG (Hill et al. 1996), each covering a period of 108 days, starting from 1995 May 7 and 
ending on 2003 August 16. As a measure of the solar activity for each data set, we use the 
mean radio flux at 10.7 cm during the time interval covered by each data set as obtained 
from the US National Geophysical Data Center (www.ngdc.noaa.gov/stp/stp.html). 

For each data set we use modes with two sets of degrees; one set consists of modes with 
degree 5 < i < 25 and we refer to this as the low-degree set, the second set consists of all 
modes with £ < 60 that have their lower turning point, > 0.715 and we refer to this 
as the intermediate-degree set. The lowest degree modes that satisfles this criterion for our 
frequency range of 2-3.5mHz is i = 30. Very few modes below i of 42 satisfy this criterion. 
The first set of modes sample both the CZ base and the He II ionization zone and hence we 
fit the entire form given by Eq. (2), while modes in the second set sample only the He II 
ionization zone and hence we fit the mode only to the first two terms of Eq. (2). We do not 
use modes with £ < 5 because they have larger errors. We do not use modes with degree 
higher than £ — 60 because the degree-dependence of the signal becomes difficult to model 
and fit. 
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The quantity we are most interested in is the amplitude of the oscillatory signal from He 
II ionization zone as a function of the magnetic activity level of the Sun when the data were 
obtained. Since the amplitude is both degree-dependent as well as frequency-dependent, we 
use the averaged amplitude in the frequency range 2 to 3.5mHz after the degree dependence 
is removed. This is the same as what was done by Basu et al. (1994), Basu & Antia (1994) 
and Basu (1997). The error on the result is determined by Monte-Carlo simulations. The 
fit to the signal from the He II ionization zone, after removal of the degree dependence, is 
shown in Fig. 2 for one set of low-degree GONG data. 

3. Results 

The amplitudes of the oscillatory signal that arises from the He 11 ionization zone arc 
plotted as a function of the solar activity index in Fig. 3. The 10.7 cm flux is in units of the 
Solar Flux Unit (SFU), i.e., lO^^^ j g-i ^^-i ^laNe plotted the four sets of data 
separately. 

We find that in all cases, the amplitude decreases with increasing solar activity, a result 
that is expected if the active-region results can be applied to the global Sun. We can 
fit straight lines to the data quite easily — the large errors do not justify fitting more 
complicated trends. We can see that for all 4 sets, the straight fine has a finite slope, however, 
the slope is not always statistically significant. This is particularly true for the GONG low- 
degree data. The MDI low-dcgrcc sets show only a marginally significant decrease. However, 
the intermediate-degree sets for both GONG and MDI show a reasonably significant trend 
(fa 4a). The scatter in the plots is somewhat less than what the errors on the points should 
suggest, so it is likely that the errors have been overestimated and significance of the slope 
underestimated. 

It is not completely surprising that the low- and intermediate-degree sets give us some- 
what different results. For the low-dcgrcc sets, wc need to correctly remove the signal from 
the CZ base to get the correct amplitude from the He II ionization zone. Simulations per- 
formed with different solar models show that this process leads to substantial systematic 
errors in the results. The intermediate- degree modes are not affected by the CZ base at all, 
and hence the signal due to the He II zone is cleaner and easier to measure. Thus we put 
more weight on the results obtained from the intermediate degree modes. 

Figure 4 shows both GONG and MDI intermediate-degree results plotted as a function 
of the 10.7 cm fiux. It is clear that both GONG and MDI show very similar trends, which is 
encouraging since these are independent projects. A straight line fit to all the points shows 
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that the slope is 5.5(7 from zero, a result that is reasonably statistically significant. We 
therefore, can conclude with a degree of confidence that the region of the He II ionization 
zone changes with change in solar activity. In particular, the magnitude of the dip in Fi, 
which is what causes the signal we are looking at, decreases with increasing activity. It 
should be noted however, that there should be a correlated noise component between the 
GONG and MDI data sets since they are observing the same object, hence the increase in 
the significance may be smaller than what we find on combining the results. A combination 
of all the results (GONG and MDI, low and intermediate-degree set) gives a slope significant 
at the 4.7(7 level. 



4. Discussion 

The easiest interpretation of a change in the magnitude of the depression of Fi at the He 
II ionization zone is a change in the abundance of helium. However, that interpretation does 
not apply in this case, since we are talking of cyclical changes over very short time-scales. 
The only change in helium that is expected in the CZ is a monotonic decrease due to the 
gravitational setthng of helium that takes place over very long time-scales. We need to look 
at the the equation of state (EOS) to understand the changes. The presence of magnetic 
fields can change the effective EOS because of contributions of the magnetic fields to energy 
and pressure. A change in the EOS changes the shape and magnitude of the dip in Fi at 
the He II ionization zone for the same helium abundance. 

Figure 5 shows the relative difference in Fi between two solar envelope models, one 
constructed with the so-called MHD equation of state (Hummer & Mihalas 1988; Mihalas, 
Dappen & Hummer 1988; Dappen ct al. 1988), and the second with the OPAL equation 
of state (Rogers, Swenson & Iglesias 1996). Both models were constructed with the same 
opacities, have identical helium and heavy element abundances in the CZ (0.242 and 0.018 
respectively) and were constructed to have the same CZ depth (O.287R0). We can see that 
the models have substantial differences in Fi in the region of the He II ionization zone and 
higher. The amplitude of the He II signal for the MHD model is 1.164/iHz and that of the 
OPAL model is 1.027/iHz. 

The difference between the He II amphtudes of the two models (0.137//Hz) is larger than 
the total range of change in amplitude seen in Fig. 4 (which is only 0.099//Hz according to 
the linear fits to the results). Thus as a first approximation we can say that the change in 
Fi near the solar helium ionization zone between solar minimum and maximum is less than 
that in Fig. 5, i.e., less than about 4%. However, one must be careful about this number; the 
Fi differences in Fig. 5 have some fairly sharp features that could be the cause of the large 
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difference in amplitude for the two models. It is unlikely that a small change in the EOS due 

to magnetic fields would cause spiky changes, and hence for 0.099/iHz could imply a larger 
difference in Fi. Inversions could easily detect differences in Fi less than the estimated Fi 
change between solar minimum and solar maximum. However, to do so reliably, and with 
any degree of statistical significance, we need mode sets with reliable and precise high-degree 
modes (Di Mauro et al. 2002; Rabello-Soares et al. 2000). 

As far as the parameter is concerned, we cannot draw any conclusions since as can 
be seen in Fig. 6, the points are widely scattered. 

One expects solar cycle related changes to have strong latitudinal dependences. We 
repeated the study described above for different latitudes, but the errors are too large to be 
able to observe any latitudinal dependence. 

5. Conclusions 

We find evidence that suggests that solar structure changes with change in solar activity 
in the layers around the He II ionization zone (i.e., 0.98 R© and thereabouts). The depression 
in the adiabatic index Fi in the He II ionization zone decreases with increasing solar activity. 
This is the first evidence to suggest that solar structure changes with solar activity. 

This work utilizes data obtained by the Global Oscillation Network Group project, 
managed by the National Solar Observatory which is operated by AURA, Inc. under a co- 
operative agreement with the NSF. This work also utilizes data from the Solar Oscillations 
Investigation/ Michelson Doppler Imager (SOI/MDI) on the Solar and Heliospheric Obser- 
vatory (SOHO). SOHO is a project of international cooperation between ESA and NASA. 
The MDI project is supported by NASA contract NAG5-13261 to Stanford University. This 
work was supported by NSF grants ATM 0206130 and ATM 0348837 to SB 

REFERENCES 

Basu, S. 1997, MNRAS 288, 572 

Basu, S. 2002, in From Solar Min to Max: Half a Solar Cycle with SOHO, Proc. SOHO 11 
Symposium, ed. A. Wilson, ESA SP-508, 7 

Basu, S., Antia, H. M. 1994, MNRAS 269, 1137 



- 7- 



Basu, S., Antia, H. M. 2000, ApJ, 541, 442 

Basu, S., Antia, H. M., ApJ, 2003, 585, 553 

Basu, S., Antia, H. M., Narasimha, D. 1994, MNRAS 267,209 

Basu, S., Antia, H. M., Bogart, R. S. 2004, ApJ, 610. 1157 

Dappen W., Mihalas D., Hummer D. G., Mihalas B. W. 1988, ApJ 332, 261 

Di Mauro, M. R, Christensen-Dalsgaard, J., Rabello- Scares, M. C. and Basu, S. 2002, A&A 
384, 666 

Eff-Darwich, A., Korzennik, S. G., Jimnez-Reyes, S. J., Prez Hernndez 2002, ApJ, 580, 574 

Elsworth, Y., Howe, R., Isaak, G. R., McLeod, C. P., New, R. 1990, Nature, 345, 322 

Gough, D. O. 1990, in Lecture Notes in Pliysics, eds., Y. Osaki, H. Sliibahashi, (Springer: 
Berlin), vol. 367, p 283. 

Hill F., et al. 1996, Science 272, 1292 

Howe, R., Komm, R., Hill, F. 1999, ApJ, 524, 1084 

Howe, R., Christensen-Dalsgaard, J., Hill, F., Komm, R. W., Larsen, R. M., Schou, J., 
Thompson, M. J., Toomre, J. 2000, ApJ, 533, L163 

Hummer D. G., Mihalas D. 1988, ApJ 331, 794 

Libbrecht, K. G., Woodard, M. F. 1990, Nature, 345, 779 

Miglio, A., Christensen-Dalsgaard, J., di Mauro, M. P., Monteiro, M. J. P. F. G., Thompson, 
M. J. 2003, in Asteroseismology Across the HR Diagram, eds. M.J. Thompson, M.S. 
Cunha, M.J.P.F.G Monteiro (Kluwer:Dordrecht) p. 537 

Mihalas D., Dappen W., Hummer D. G. 1988, ApJ 331, 815 

Monteiro, M. J. P. F. G., Christensen-Dalsgaard, J., Thompson, M. J. 1994, A&A, 283, 247 

Monteiro, M. J. P. F. G., Christensen-Dalsgaard, J., Schou, J., Thompson, M. J. 2000, in 
Helio- and asteroseismology at the dawn of the millennium, Proc. SOHO 10/GONG 
2000 Workshop, ed. A. Wilson, ESA SP-464, 535 

Kosovichev, A. G., Duvall, T. L., Jr., Scherrer, P. H. 2000, Solar Phys., 192, 159 



-8- 

Kosovichev, A. G., Duvall, T. L., Jr., Birch, A. C, Gizon, L., Scherrer, P. H., Zhao, J. 

2001, in Hclio- and Asteroseismology at the Dawn of the Millennium, Proc. SOHO 
10/GONG 2000 Workshop, ed., A. Wilson, ESA SP-464, 701 

Rabello-Soares, M. C., Basu, S., Christensen-Dalsgaard, J., Di Mauro, M. P. 2000, Solar 
Phys. 193, 345 

Rogers, F. J., Swenson, F. J., Iglesias, C. A. 1996, ApJ 456, 902 
Schou J., 1999, ApJ 523, L181 

Vorontsov, S. V. 2001, in Helio- and Astero-seismology at the dawn of the millennium: Proc. 
SOHO 10/GONG 2000 Workshop, ed. A. Wilson, ESA SP-464, 563 

Vorontsov, S. V., Christensen-Dalsgaard, J., Schou, J., Strakhov, V. N., Thompson, M. J. 

2002, Science, 296. 101 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-9- 
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Fig. 1. — The adiabatic index Fi for a standard solar model and two active regions. The 
region AR9893 is much stronger than AR 8040 and their magnetic indices as defined in Basu 
et al. (2004) are shown in the figures. The Fi profiles for the active regions were determined 
from results of Basu ct al. (2004) assuming that the Fi profile of the quiet regions of the Sun 
is like that of a standard solar model. The results would not change appreciably if, instead 
of a standard solar model, Fi obtained from inversions of quiet-sun data had been used. 
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Fig. 2. — The fourth-difference of frequencies of a set of GONG data with the degree- 
dependence in the signal removed. Also shown is a fit to the signal in the Fi depression. 
Data for frequencies with 5 < £ < 25 are shown. 
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Fig. 3. — The amplitude of the oscillatory signal averaged over the frequency range 2-3.5 
mHz plotted as a function of the 10.7 cm radio flux, which is an indicator of solar activity. 
The 10.7 cm flux is plotted in units of 10~^^ J s~^ m'^ Hz— 1 (Solar flux units). The different 
panels show the results of the different degree ranges of the GONG and MDI sets. The linear 
least squares fit to the points is shown as a continuous line. Each panel also lists how many 
standard deviations away the slope of the fitted line is away from zero. 
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Fig. 4.— The GONG and MDI results obtained for modes with ^ < 60 and r* > O.TISRq 
plotted as a function of the 10.7 cm flux. The black continuous line is a least-squares fit 
to all the data points. The dotted lines show the la spread because of errors in the fitted 
parameters. 
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Fig. 5. — The relative difFerence in the adiabatic index Fi between two solar envelope models, 
one constructed with the MHD and the other with the OPAL equation of state. The models 
were constructed with the same helium and heavy metal abundances and the same CZ depth. 
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Fig. 6.- 



The variation of as a function of the 10.7cm flux. 



